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ABSTEAUT 


Incoherent infrared scattering of silicon is reviewed and scat- 
tering from a crystal grown along a <110> axis has been examined. It 
is concluded that the incoherent infrared scattering arises from regions 
in the cyrstal bounded by dislocation loops. The observed radii of 
cyration of the scattering regions, +5 As » agree with the size of 
dislocation loops which are observed using Dash's decorating technique. 
Analysis of the small angle scattering precludes the presence of oxide 
particles smaller than about 1M, and microscopic examination discloses 
no particles larger than this indicating that there has been no coalesence 


of oxides. 
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CHAPTER I 
INTRODUCTION 


Infrared radiation has been employed by many investigators in at- 
tempts to study imperfections in the structure of silicon crystals. 

The majority of efforts have been directed toward the deliniation of 
dislocation and impurity defects. Chemical etching has long been used 
as a means of visual detection of the emergence of dislocations at 
crystal surfaces. A method of creating images of dislocations within 
the bulk of the crystal by a decorating technique has been developed by 
Tees at the General Electric Laboratories. He was able to trace 
visually the dislocation lines through the crystal and detect their e- 
mergence at etch pits at the surface. X-ray extinction microscopy has 
been applied in determining dislocation NenancHeen oe A unique method 
using an integrating sphere to detect light scattering by both disloca- 
tions and oxide impurities was used at the General Telephone and Elec- 
tronics Laboratories. 

This investigation was directed toward the application of incoherent 
scattering techniques as a means of Nalseielg size determination in sili- 
con crystals. Small angle X-ray diffraction is a well established 
method used to measure particle sizes of the order of magnitude of X-ray 
wavelength, However, this method is not applicable to structure deter- 
minations of sizes large compared to the wavelength of the radiation used. 
Theoretically, with longer wavelengths, one could determine the sizes of 
larger particles. The problem encountered is that softer X-rays and 


visible light are attenuated to quickly by matter. The wavelengths in the 


near infrared spectrum are of the magnitude of dislocation and impurity 
networks observed by the methods referred to above. Therefore, in- 

coherent scattering of infrared appears to offer possibilities as an- 

other means of determining the dislocation densities within silicon 
crystals. It is also proposed that a diffraction system can be built easily 
and cheaply thus avoiding the highly designed and costly equipment associ- 


ated with X-ray diffraction and small angle cameras. 





Cire “Ti 


REVIEW OF DISLOCATION AND OXIDE IMPURITY DETECTION IN 


STLICON CRYSTALS 


1. Chemical Etching - A method of identifying dislocations in LiF 
crystals by observations of etch-pits reported by euuiinare has been 
applied to semi-conductor materials by other investigators (Vogel et ah 
renee). Dash's work was devoted to the study of dislocations in sili- 
con crystals. He reported etch-pits several hundred microns in diameter 
extending a comparable distance into the crystal when etched overnight 


with an echant consisting of the following parts by volume: 1 50% HF, 


3 70% HNO 10 glacial acetic acid. One can obtain the dislocation 


3? 
density of a crystal by counting the etch-pits which mark the emergence 
of dislocations at the crystal surface. Dash observed deep pits on any 
surface independent of its crystallographic direction. The disloca- 

tion density within any one crystal may vary strongly with orientation, 
and etch-pit counts on only one surface can lead to erroneous conclusions 
about crystal perfection. sehwuceker reported errors up to 100 per cent 
in crystals of low dislocation density. Etch-pit counting is universally 


applied by investigators as a means of obtaining dislocation densities in 


silicon crystals, 


2. Decoration - Me oS describes a technique for observing visually 
the dislocation lineage within the bulk of silicon crystals. Copper 
atoms diffused into the crystal were precipitated in aggregates along 
dislocations after appropriate heat treatment. Silicon but not copper, 


is transparent to neai infrared radiation. Using an infrared image tube 
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fitted with a microscope he was able to obtain images of dislocation 
lineages in the crystals. A one to one correlation between the images 
and dislocations was obtained by tracing the lines to etch-pits at 
crystal surfaces. Precipitates are frequently not found on pure screw 
dislocations, but when formed they appear as tightly wound helicies. 

In <10> grown crystals he observed dislocation loops ranging up to 

50 microns in diameter. Copper atoms appear as a fine, speckeled sub- 
structure in dislocation-free crystals. This is apparently a result of 


agglomeration of the precipatate throughout the crystal body. 


2. Light Scattering 


a. Integrating Sphere Technique - Sinema reports a method of 
detecting incoherent scattering of infrared by "bad regions" within 
silicon crystals. He calculated the number of scattering centers (re- 
gions outlined by dislocations) necessary to produce detectible scatter- 
ing. Three assumptions were made: 1. Imperfection scattering, to be 
observable, must be ten times greater than thermal scattering,. 2. The 
scattering centers are spherical in shape, 3. The radius of the scatter- 
ing centers is of the order of 0.1 micron. Since the particle size is 
small comparted to the wavelength of the radiation ( = 1.1 ), ex- 
pressions for Rayleigh scattering are developed relating intensity to 
the size and number of bad regions per unit volume. Theoretical consider- 
ations suggest that a dislocation density of at least 10° per mS is 
necessary to produce observable scattering. This was confirmed in 
Schwuttke's experiment, but the results do not constitute a confirma- 
tion of the assumed particle ae The object of the study was to 


determine the dislocation density necessary to produce scattering. 
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As previously noted, Dash reported dislocation loops of from 2-50 micrcens 
in diameter. Also the loops tend to lie in favorably orientated planes 
thus suggesting a platelet configuration. Regions of 0.1 micron (the 
assumed particle size in this study) described by dislocations cannot be 
resolved by Dash's method, but since the scattering depends on a radius 
of gyration there exists an infinite number of particle size and shape 
combinations which could produce the same scattering effect. Lacking 
Supplementary information, one cannot draw definite conclusions about 
either size or shape. 

Schwuttke's ingenious experimental technique of observing very small 
intensities of scattered radiation warrants a brief description. The 
optical arrangement of the instrument appears in figure 2-1, 

A gold-plated integrating sphere contains a tungsten filiment light 
source S. The black body B and openings 0, and 0. define a conical umbra 


into which no light originating inside the sphere can penetrate. How- 


ever, scattering by a crystal placed outside of 0, will permit light to 


Zz 
enter the "dark'' area. With a suitable optical system and infrared image 


tube, scattered light can be observed. 


b, X-ray Extinction - X-ray diffraction techniques have been used 
to substantiate dislocation counts in crystals obtained by the previous- 
ly described techniques, The double-crystal spectrometer rocking curve 
method is commonly applied to determine crystal perfection. fang re- 
ports this method to be unreliable for dislocation densities below about 
10° per aie X-ray topographs can present an unambiguous display of 
dislocation lineage within a crystal. He was able to observe directly 
individual dislocations by copper atom decoration. Lang's results sub- 


stantiate the preference of the {119> direction for straight line 
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Fig, 2-1 Schwuttke's Scattering Integrating Sphere 
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dislocations teported by Dash. Arrays of hexagonal loops are formed 
by the closing of <1106> edge dislocation lines by pure screw components, 
This suggests a thin two dimensional platelet type of configuration for 


a particle shape, oe 


c. Additional Results - One ircomplete report of dislocation 
observation was found in the literature. The results, but apparently 
not the experimental procedures, were published. Patel, Wagner and occ 
have observed the aggregation and collapsing of supersaturated vacancies 
to form condensed prismatic dislocation loops upon cooling silicon crystals 
from the melt. They report that indirect evidence of this was obtained 
by Chandhuri and Patel from plastic deformation experiments. The dis- 


location loops observed after ceformation were estimated to be about 1 


micron in diameter, 


4. Oxygen Impurities - Oxygen impurities have been observed by both 
decoration and light scattering techniques. Several investigators have 
reported the effect of oxygen in contributing to erroneous dislocation 
etch-pit counts but have not used echants as a means of oxide detection. 
Debye and Megermandier = in an incomplete report state they have 
obtained Si0, particle sizes by measuring angular dependence of scatter- 
ed infrared radiation. They found particle sizes of the order of the 
wavelength of the light. The sizes are strongly dependent upon the heat 
treatment given the specimens. eer reports the effect of long (50 
hours) heat treatment at 1000°C on the optical properties of silicon. 
a to heating, silicon exhibits a sharp infrared absorption band at 


a wavelength of 9 microns. The absorption is proportional to oxygen con- 


centration. Following heat trcatment, the band decreases but without an 
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accompanying lowering iu oxygen concentration. He attributes this to 
the clustering of silicon oxide into particles of lower dielectric 
constant than the silicon matrix. 

Kaiser also makes an attempt at S10, particle size determination 
by observation of Rayleigh scattering. He assumes a spherical particle 
of size smaller than the wavelength of radiation employed. Since Si and 
Si0, do not absorb at wavelengths of 1-3 , changes in intensity will 
be due to Rayleigh scattering and not absorption. For spherical parti- 


cles, the absorption cvefficient (CC) due to Rayleigh scattering is 


oo = KN, V7\4 Dey 


where K is a constant, Np the number of scattering centers, V the total 
volume of scattering centers and A the wavelength. Following the in- 
verse 4th power law, a log-log plot of A vs. A should yield a straight 
line of slope equal to -4. Experimental results produced a slope of -3.4 
which adequately substantiates the existance of Rayleigh scattering. 


From 2-1 and the auxiliary equation 


No Mo 
Y= —_— ye yey) 
Ne ( 288 


where Ny is the total number of oxygen atoms per aa M the molecular 
wright of Si0,, Q the density, and A avogradro's number, the number of 
Scattering centers, Np can be calculated. The number of oxygen atoms 
per cluster, No/Np, is obtained and from the assumption of spherical 
particles he calculates the diameter of scattering to be O.1 Mr. This 
value is consistant with the magnitudes observed by Debye and Lederhand- 


ler from angular dependence scattering experiments. 





Schwuttke’s vork produces the most conclusive and significant 
information about the distribution of oxygen in silicon. He calculated 
that a minimum of fo" oxygen atoms per “ii was necessary to produce 
scattering. This was confirmed by observations with his integrating 
sphere technique. An investigation was then conducted to determine 
the distribution of the oxide in the crystal. By application of etch-pit 
counts, decoration and light scattering techniques to the same specimens 
he was able to make observations of the effect of dislocations upon the 
oxide particles, It was found that in decoration by copper atoms, the 
precipitation is influenced by the oxygen content of the crystal. This 
suggests that oxygen is present in the dislocation networks and acts as 
nucleation centers for the copper precipitates. The copper prefers re- 
gions of high oxygen content thus delineating high density oxide clusters. 
Dislocations offer favorable sites for Si0, clusters. In relatively dis- 
location free crystals, clustering by heat treatment tends to take place 


in the bulk of the crystal. 


>. Summary - From the foregoing review the following conclusions are 
presented: 

1, Dislocation networds will produce incoherent infrared scatter- 
ing and allow one to make qualitative conclusions about their distribu- 
tion and size. 

2, The assumption of a spherical shaped particle size for dis- 
location regions is not necessarily valid. Decorating techniques suggest 
a fishnet structure with a particle shape aprroaching thin platelets. 

3, Particles formed by dislocation networks have some degree of 


preferred orientation. 





4, Oxygen tends to concentrate on dislocation sites. leat 
treatment at lOOO°C produces higher concentrations of Sid, on these 
sites through clustering. This constitutes an increase in particle size. 

5. In crystals of low dislocation density, Sid, molecules will 
cluster in the body of the crystal. Assumptions of a spherical cluster 
shape are not substantiated. 

6. Dislocations and oxide impurities of appropriate size ere 
present to cause incoherent scattering. 

7. Direct evidence exists of dislocation loops as large as 20 A 
in viareter. Strong indications exist for sizes as small as La dia- 


meter, 
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CHAPTER IIT 
THEORY OF SMALL™ANGLE SCATTERING OF X-RAYS 


A. General - The purpose of this study is to investigate the feasi- 
bility of applying small-angle scattering techniques to the determina- 
tion of dislocation networks or impurity concentrations in silicon 
crystals, These heterogeneities constitute a particle distribution 
within the material, Therefore, it is appropriate that a review of the 
theory of small-angle scattering by particles be presented at this point. 
The early theory was primarily developed on a qualitative basis by 


13,14 and later extended to yield more quantitative results by 


Guinier 
Guiner and Fouenee Small angle scattering is a phenomenon commonly 
observed with visible light. The halo observed around the moon due to 
Scattering of the light by small mist droplets is the result of small 
angle scattering. To observe such an effect, the size of the scattering 
particles must be of the order of ten to a hundred times the size of the 


wavelength of the scattered light. For this type of scattering the 


following approximate relationship is observed: 


c = oT (3-1) 


in which € is the angle of diffraction, X the wavelength, and d the 
mean diameter of the particles. The similarity between this relation 


and Bragg's Law, 
ee 2d Sin S (3-2) 


is apparent. However, one must not associate this type of diffraction 


11 





with that produced by the regular periodic cyrstallites to which Bragg's 
Law is applied. It is, in fact, the very nature of the variation, or 
non-periodicity of structure which gives rise to the small angle scatter- 
ing to which the theory will be applied. Guinier extended the relation- 
ship in equat ion (2-1) to include the same type of diffraction phenomenon 
in the region of X-rays. If the wavelength and particle diameter were 
both reduced in the same proportion, particles of the order of ten to a 
hundred angstroms should likewise produce detectable small angle scatter- 
ing. 

X-rays are scattered by the electrons of the diffracting medium. 
Ordinary X-ray patterns are obtained from the scattering of electrons 
of atoms situated at regular coherent lattice points. However, low- 
angle scattering depends only upon the external dimensions of the parti- 
cle, that is, the contour of the change in election density and not its 
internal structure. The X-ray diffraction of this type ordinarily ob- 
served is of the form of continuous scattering at angle of two degrees 


and below. 


B. Theory - In order to make application of the phenomena to the deter- 
mining of particle size, a relationship between scattered intensity and 
angle containing some dimensional parameter characteristic of the parti- 
cle must be developed. This parameter obviously must be sensitive to 
the difference in the electron density due to the inhomodenity created 
by the particle, The derivation of appropriate formulas will be found 
in Appendix I, and only the results presented below. The characteristic 
particle dimension is designated by the symbol R. It must be dependent 
upon the electron distribution which in turn is a function of the atom 


distribution within the particle. The total number of scattering points 


is therefore N atoms times the atomic number. The relative displacenent 
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of these points from sorie reference, in particular the center of gravity, 
is of interest. Take R? as the mean of the squares of the distances of 
each atom from the center and attach to it a value equal to its atomic 
weight. The value of R then is exactly that quantity in mechanics call- 
ed the radius of gyration. It is a distance such that MR? = I where M is 
the particle mass and I its moment of inertia relative to the center of 
gravity. This parameter therefore is only an indication of mass distri- 
bution and not shape. It is because of this very point that the basic 
theory still leads to only qualitative or semi-quantitative results, 
The majority of Guinier's later user was devoted to the refinement of 
this aspect. 

Appendix T yields the relationships: 


for orientated particles ...... 
> exe(- 41” pic? 
t/t. Sealy ECuxe ae (3-3) 
for randomly orientated particles ...... 
+3 
= z . 41T 272 (3-4) 
t/t, = Nn exp =a RE } 


in which T = scattered intensity over an angle € (in radians) 


I = the scattering intensity by one electron 
(obtained by Thomson's formula) 


N = the number of particles 

n = the difference in the number of electrons in the 
particle and an equal volume of the external 
medium, 


Taking the logarithom of each side of(2-3)one obtains: 


2 
oe ee egies 
Logi ol = Log Nn r. = 5 1.08 1 eR & 





3 
A pidt of Ue, io: Nee & ahouwls Ste dd a stracght line for small 


Dall _ ee 2 
values of € ', its slope, w, being ~41T Z AX Loy io = ek. his 


yields a value of 


R= Atm = 6 Woes NVM 


41? Loq,,€ 





for orientated particies and .418 AN for randomly orientated 
particles. Usually one plots Log jot as a function of the square of 


the magnitude of the vector h, where: 


4ir Sin © Ay ae 


In) = i. 
A a for small angles 
(9 = 1/2 the scattering angle) 
C. Interpretation of results - A straight line curve predicted by 


equations (3-3) and (3-4) is seldomly obtained in practice. To obtain 

a strictly linear curve demands only one discrete particle size and of 
course good precise experiinental data. Figure 3-1 is a good Vxample of 
these conditions for a solution of albumin. The radius of gyration is 
obtained from the slope of the linear plot. From other information about 
the particle such as the density, modecular weight, and chemical struc- | 
ture, R can be transformed into a particle shape with specific dimensions. 
Now take the case presented in Figures 3-2 and 3-3, What statements can 
be made about particle sizes for these experimental data? The theory of 
small angle scattering exhibits its greatest validity for particles of 
dilute concentrations. This is due to the basic hypothesis that it is 

the difference in electron density of the particle and the medium upon 
which such scattering depends, For increasing concentrations and orienta- 
tions a point is reached in which there appears to be a degree of long 


range order and particle interference, This tends to create a step or 
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Fig. 3-1 Low angle scattering of a solution of albumin,after 
Guinier 
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Fig. 3-2 Particle interference Fig. 3-3 Mixtures of particle 
sizes 





even a maxima in the otherwise continuous scattering curve. Figure 

3-2 represents this situation. The most representative of scattering 
patterns is represented in Figure 3-3. Only qualitative information 

can be obtained from these curves. In curve l, the variation of the slope 
at low angles is sufficiently small to permit an extrapolation to the 
origin. The slope at the origin gives an average radius of gyration. One 
can say that the sample does not contain particles of too large an R to 
escape measurement. Curves 2 and 3 do not permit such an extrapolation. 
The conclusions to be made are: Curve 2, the sample contains many large 
particles and few small ones; Curve 3, the sample contains some very large 
particles but not many small ones. 

Guia? and meyer’ extend the theory in an attempt to interpret 
quantitatively curves exhibiting mixtures of particle sizes, His analy- 
sis is based upon distribution functions in which he applies Maxwellian 
distributions. A more straight-forward graphical procedure has been 
developed by Jellinek, Solomon and fFadtéucWen’ in which the contribu- 
tions of particle size groups are successively subtracted out of the 
curve. The curves of Figure 3-3 represent such a sum of the contribu- 
tions of many particle sizes, They point out that the scattering of X- 
rays by a homogeneous assemblage of spherical particles of radius R can 


be given by the equation, 

ee 
2 .O.221h R - oe 
T(e) = C¢wre 5) 
where [I is the scattered intensity as a function of @ (half the scat- 


4x Sin€ 
oN 


scattering and C is a constant for a given material. The distribution 


tering angle), h = , W is the weight of tle specimen doing the 


of particle sizes may be represented by the Function W(R) where W(R) is 
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the weight fraction within the sample of radius R. Summing the total 
scattered intensity leads to the expression for the curve: 
ca 2 =—GOraet }? ae 
rle)=cf W (R)R™e dk 3-6) 
a) 

The curve of Log I (8) vs, ne will, at the higher values of eo yield 
a slope representing a single value function of R and with an intercept 
equal to CWR?, Equation (2-6) suggests possibilities of obtaining W(R) 
as a continuous function of I(0) by application of the Fourier transform 
method. A treatment of this analysis is given by Becaeas involving a 
great detailed process. The Jellinek, Solomon, and Fcankuchen method 
greatly simplifies matters by obtaining in place of the continuous func- 
tion of W(R) a graphically obtained set of discrete fractions of W(R) over 
the observed range. Their data and graphical analysis for alumina gel 
are presented in Figure 3-4. The tangent to the curve at the greatest 
observed angle of scattering is taken and extrapolated back to the 


ordinate as the intercept K A point by point subtraction of the 


1° 
tangent from the curve yieldsthe new curve (dashed in the figure). This 


represents the scattering contributed by all particles of size larger 


than the value yielding slope m From this curve the tangent m, is 


1° 2 


taken and the intercept Ky is obtained. The process is repeated until 
extrapolation and point by point curve subtraction is no longer feasible. 
The slopes are calculated and from an appropriate graph of slope vs. R 
(Jellinek, Solomon and pranieecheane the discrete radii are obtained. 
Since the intercepts are propcrtional to Rw(R), the quantity K,/ Ro 


: .th ; : ae 
represents the fraction of the 1 particle contained within the sample. 


The complete calculation for carbon black is presented in Table l. 
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TABLE 1 


Calculation of Particle Sizes from Small-Angle Scattering Data for 
Alumina gel (Jellinek, Solomon, and Frankuchen/) 


Ordinate Calculation of Slope 2 
sel Ae y 2 Ke yiot Weight 


Intercept log I r log I r R R Proportion 
Ry 0.61 0,484 630 0.000769 10 1,000 6.10 0.194 
K 4,3 ja oo SN Js 7 0.00358 Ze 10,600 4.06 0.129 


K 39 2.591 240 0.0108 37 DO,600 “Wye ie 0. 245 


3 
K, 115 3.061 140 0.0219 54 157,000 7.32 0.232 
Ke 140 3,146 107 0.0294 62 238,000 5.88 0.186 
Ke 48 2.681 33. 0.0812 103 1,092,000 0.44 0.014 





Total 31.51 1,000 


Jellinek et al in their original publication acknowledged the 
apparent empirical nature and only piecewise continuous size distri- 
bution of this method of analysis. They could only defend its validity 
by stating that it works for the systems (alumina gel and carbon black) 
to which it was applied. Ina later analysis of alumina a they 
were unable to obtain consistent results. 

Bragg, Corwin and aucerey investigated mixtures of silica sols of 
two known particle diameters in an attempt to verify methods of particle 
size distribution, They were able to substantiate the prediction that 
the scattering at small angles arises predominantly from the larger 
particles and that the smaller ones control the scattering at larger 
angles. The application of Jellinek's graphical analysis was unsuccessful 
in reproducing the known particle sizes introduced, They state, however, 


that much of their scattering data was obtained at values of h beyond 
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which the Guinier approximation is valid. A method due to Porod 
which relates scattering at larger angles and total surface area 


was also applied. Porod predicts the following relationship: 
ry 
T(h) = 2z7re Y;4 (3-7) 


in which Qis the mean electron density, S the total particle surface 
area and h as previously defined. In collimating systems of large 

slit widths, the quantity is replaced by h me plow ofMitr [vs re h 
should yield a straight line of slope equal to -3 (using i) and inter- 
Caperof 2x Qs For spherical particle disperse system the geometri- 
cal relationship R= 3 V/¢ , where V is the particle volume, obvi- 
ously holds, An analogous surface averaged particle radius can be 
obtained for polydispersed systems. The calculated radius favors those 
particles which make the largest contribution to the surface area, 

The above survey substantiates Guinier's original prediction that 
the application of small angle scattering in obtaining accurate quantita- 
tive results is indeed limited. For widely separated spherical particles 
of one diameter, exact sizes can be measured, Nou-spherical particles 
can be quantitatively investigated if complementary information such as 
volume, surface area, and molecular weight are known. Mixtures of sizes 
and/or unknown shapes lead to qualitative results of widely varying 


validity depending upon the means of analysis employed. 





CHAPTER TV 


APPLICATION OF SMALLWANGLE INFRARED SCATTERING TO PARTICLE SIZE 


IN SILICON CRYSTALS 


A. General - Guinier developed his theory of small-angle scattering 
by an extension of the classical theory of atmospheric scattering of 
visible light down to the region of X-rays and particles of comparable 
size, It therefore seems reasonable that a similar extrapolation can 
be made toward the region of the infrared spectrum and particle dimen- 
sions of the order of microns. The limitation to materials which are 
relatively transparent to infrared radiation is obvious, An inspection 
of equation (3-1), ex 4 suggests two possibilities. One might at- 
tempt to increase the wavelength, maintaining a given particle size, 
and thus obtain scattering at larger angles. This would eliminate the 
complexities associated with small angle cameras and diffraction tech- 
niques. However, in the case of X-rays, longer wavelengths are increas- 
ingly attenuated by matter. This introduces problems associated with 
reduced intensities. The theoretical limitation placed upon this first 
alternative is that the expressions developed in the small-angle theory 


are derived using an approximation called the “law of Guinier". Equa- 
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tion (3-3) and (3-4) for the scattered intensity contain the term e 


in which K represents =e . As shown in Appendix I this exponential is 


an approximation of the series 
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The approximation coincides with the exact expression only up to the 
4th power of h. Therefore, for large h, the relationship 

t/z. = shale ©. he R cs 
is no longer valid. In fact, the Guinier analysis breaks down for hR &1. 
The above approximation is the law of Guinier, Therefore, one cannot 
increase the wavelength in an attempt to obtain scattering at larger 
angles and still apply Guiner's small-angle scattering theory. 

The other suggested possibility as previously mentioned is that of 
increasing both X and d, maintaining the same proportion. In this 
manner, the angle limitation placed upon € is not exceeded and larger 
size particles can be measured. The study of silicon heterogeneities 
will be directed toward using this relationship. 

Kaiser, Keck and vanes were able to show by means of a 9,1 micron 
wavelength absorption band the existance of Sid, clustering in silicon 
crystals, In Chapter II it was pointed out that Schwuttke was able to 
make some semi-quantitative conclusions of oxygen impurity clustering 
and dislocation networks in silicon by infrared scattering. Debye and 
feerhandier-° stated that they were able to obtain information on Sid. 
particle size in silicon by measuring the angular distribution of infrared 
scattering. None of the above investigators were able to establish a 
discrete particle size. 

The problem of interpretation of scattering data still remains open. 
With what kind of particle shape can one associate dislocations? Can 
Si0, clustering be considered to consist of spheroids, elipsoids or 


2 


platelets? There is strong indication that oxygen precipitation occurs 
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along dislocation networks outlining a stringy network, or possibly, 


along rods of an average diameter but "infinite" length. 


B. Dislocations - Dash's decorating techniques allows one to draw some 
conclusions about dislocation networks. Grown-in dislocations tend to 
form in favorable orientations. The formation of a hexagonal structure 
in the (111) plane by the intersection of edge dislocations and screw 
components was pointed out in Chapter II. Plastic deformation induces 
dislocation loops from Frank~Reed sources in the (111) plane. Networks 
of loops ranging from 2 to 50 microns in diameter were detected. These 
observations suggest the use of equation (3-3) for orientated particles 
for low angle scattering analysis. The fulfillment of the widely separ- 
ated particle criterion can be assured by an inspection of the scattering 
profile. A dense system will exhibit particle interference and a drop 
off in intensity at zero angle. A mixture of sizes will be indicated by 
the concave upward nature of the plot of Log I vs € es 

The extension of small~angle scattering theory to the measuring of 
dislocations networks offers one gross possible shortcoming. The scat- 
tering is due to the difference in the electron density of the particle 
and the surrounding medium. Does a dislocation meet the description of 
a particle in this aspect? Dislocations by definition do represent dis- 
continuities in atomic structure. This represents only a variation in 
electronic density covering a magnitude of only atomic distances. How- 
ever, there is a stress field surrounding the core which falls off ac- 
cording to wee r being the distance from the core. Compared to the 


magnitudes of atomic forces, this may be considered as a relatively 
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long range force. From the cunsideration of Hooke's Law, this stress 
is only a manifestation of atomic displacements. aise indicates 
that one-half of the total strain energy surrounding a dislocatior lies 
within a region near the core, This offers possibilities of contribut- 
ing sufficient variation in the electron density to affect scattering 
of the infrared radiation. Furthermore, Guinier emphasized that ''the 
low-angle scattering depends on the external dimension of the particle 
and not on its internal structure."’ One can thus draw the conclusion 
that a network of dislocations warrants the description of a group of 
particles, 

As reported in Chapter II, Schwuttke has reported the detection 
of infrared scattering in silicon crystals containing dislocation 


2 4 8 Z e e 
densities of at least 10 per cm. This observation favorably supports 


the "particle nature" of dislocation networks suggested above. 


C. Oxygen Impurities - As previously noted in Chapter II, oxygen ex- 
ists as clusters in silicon in the form of Si0,. Schwuttke reported 

; ; ee 17 3 
scattering in crystals containing at least 10 oxygen atoms per cm, 

: : : 10 

Kaiser reported Rayleigh scattering at concentrations of 5X10 scatter- 
: 18 3 ; 
ing centers (10 oxygen atoms) per cm, Both investigators made the 
assumption that the si0, clusters are spherical in shape. This assump- 
tion has not been verified. A great deal of evidence leads one to be- 
lieve that precipitation of the oxide occurs at dislocation sites as 
do the copper atoms observed in Dash's decorating technique. Micro- 
scopic examination fails to reveal the presence of a separate S10, 


phase. This however only serves to insure that clusters greater than 


approximately 0.1-1 micron do exist. Submicroscopic spheriodal particles 
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and dislocation site precipitates could still be formed. The exis- 
tance of submicroscopic spheriodal clusters would obviously be most 
favorable to the small angle scattering method. They would tend to be 
nearly uniform in size and widely separated in distribution. All three 
of these conditions (spherical, uniform, dilute) characterize the 
optimum conditions for the application of the Guinier analysis. 

The alternative of the diffusion of oxygen to dislocation sites 
still exists. It is known that dislocation sites play a significant role 
in the nucleation of Excite Ceeeee The problem reduces to: once 
nucleation has commenced, what shape do the particles take? Continuous 
precipitation along dislocation lines establishes the same question of 
meeting the particle criterion as discussed in section B. 

The foregoing discussion presents strong evidence of the applic- 
ability of small-angle infrared scattering techniques to silicon crystals. 
With the appropriate complimentary information, a good estimate of parti- 
cle size and distribution should be obtained, Chapter VI will be devoted 


to the analysis of experimental data by use of the ideas presented here. 
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CHAPTER V 
BXPERTMENTAL PROCEDURE 


A. Specimen preparation - A 1,5 cm, diameter silicon single rystal 

was obtained from the Shockley Transistor Laboratory. This crystal 

was grown in the <1009 direction and measured 10 ohm-cm_ resistivity. 
Several slices ranging from 1.5 to 6 mm. were cut perpendicular to the 
growth axis and polished by standard metallographic procedures to minimize 
surface scratches and pitting. A check of the finish at 400 X magni- 
fication revealed a relatively small density of surface pitting. 

It was initially anticipated that absolutely parallel plane sur- 
faces would be necessary in order to obtain consistant results. Re- 
peated attempts were unsuccessful in cutting parallel faces. It was 
decided that a small amount of wedge could be tolerated if its gradient 
is parallel to the plane of detector rotation, If the taper were not in 
the plane of rotation, the refracted beam would contain a vertical com- 
ponent which would vary across the face of the radiation detector, A 
taper in the plane of rotation merely caused a refraction which estab- 
lished the non-diffracted or "zero angle" beam, I: at a point other 
than directly in line with the incident rays. Accordingly, the crystals 
were mounted and orientated by use of an optical gonimeter, The mea- 
sured deviation from parallel averaged about 1 degree. This amount of 
wedge was checked by means of the infrared beam using the relationship 


for minimum deviation of refracted light. For minimum deviation: 


Sin /2 (x Bs sin) 
Sin yy oh 


26 





n = index of refraction of the Silicon tor the 
wavelength used (3,449 for A = 2.3 ) 26 
A= angle between crystal faces. 


Smin = deviation angle, 
The refracted angle, S min, as measured by peaking intensity, agreed 


to within 10’ arc of the value determined by the optical goniometer. 


B. Optics - A photograph of the experimental set-up appears in figure 
5-1. The basic instrument is an American Optical Company table spectro- 
meter which had previously been modified for use as an optical gonimeter. 
ta an instrument used to measure small angle scattering there is a pro- 
blem of achieving adequate intensity without excessively sacrificing re- 
solution. One must make a decision as to the poorest acceptable resolu- 
tion and then strive to obtain sufficient intensity consistant with the 
sensitivity of the detector and radiation source. Two pairs of slits were 
used; one pair to collimate the incident beam and one pair to shield the 
detector from all radiation other than the transmitted diffracted beam, 
Slit separation was limited to 13 cm. by the collimator barrels of the 
spectrometer, Initially, the source of radiation used was a carbon arc. 
The carbon are is a rich source of near infrared, but difficulty was en- 
countered in maintaining a stable, well defined vapor crater which 
created alignment problems. A zirconium arc was then substituted for 
the carbon arc with a moderate decrease in intensity but considerable 
gain in stability. Several runs were made at varying slit openings and 
an acceptable signal was obtained with a width of 0.79 mn. 

The wavelength employed was governed by the optical transmission of 


silicon and the type of detectors available. Figures 5-2 and 5-3 represent 
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the factors governing the wavelergii: selected, A wave length of 2.3 mi- 
crons was decided upon because of the high peaking sensitivity of PbS 
photo-conductors at this wavelength. This choice is consistant with 
the broad transmition band of silicon in the near infrared region. The 
beam was monochromated with a “spike'' interference filter. As previous- 
ly mentioned, obtaining sufficient intensity in the beam was a problem. 
Figure 5-4 shows how the total diffracted energy is successively dimin- 
ished from that of the initial source, The energy radiated is represent- 
ed by the area under the curve of I=F (A). A detector responds to 
total energy received within the range of its sensitivity. The filter 
greatly diminishes the energy received by the detector. Additional re- 
ductions in energy are caused by components in the optical train. The 
The filter transmits only approximately 80% of the incident radiation 
at the center of the wave length. Another fraction of this transmitted 
beam is absorbed by the crystal. 

An attempt was made to increase the light intensity by concentrating 
the beam with a lens system. It was suspected that the added attenuation 
due to the glass oe ical elements would be more than offset by the beam 
concentration. Also, if the system were successful, NaCl or KBr lenses 
could be substituted for glass, thus reducing the attenuation. The 
focusing of the light with lenses did increase the detector response. 
However, this gain was obtained at the Weense of collimation and associ- 
ated directional resolution. One cannot increase the energy radiated per 
area and still maintain parallel rays. In brief, the parallel rays can- 
not be squeezed closer together to increase the photon density. To do so 


would appear to be contrary to the second law of thermodynamics. A lens 


30 


—_—-«- 
=> © > 
oa al © ound] 


—_— ewes oe ' 





system would serve only in moving the position of the source as a matter 
of convenience, The source has a finite size and cannot be concent rat- 
ed, Therefore, an increased signal was obtained through increased slit 
width and moving the source as close as practicable, consistant with 


adequate resolution. 


C. Electronics - Associated with the problem of the low intensity was 
the problem of signal amplification. Figure 5-5 is a schematic sketch 

of a typical detector circuit. In photoconducting cells, it is necessary 
to pass a direct current through them in order to record the change in 
resistance due to photo-conductivity, Maximum sensitivity is obtained in 
PbS cells with an applied voltage of the order of 100 volts. "Dark" 
resistance of a photo-conductor is usually 1 megohm. For a maximum 
transfer of power, the output or load resistor Re should be matched to 
the detector resistance and thus also 1 megohm. This creates a voltage 
divider providing approximately 50 volts across the load resistor. With 
the small signals obtained, of the order of millivolts, it was not possible 
to separate the variation in voltage from the voltage measured. Increas- 
ing the load resistor or using d.c. amplifiers merely increased the V a- 
long with V maintaining the same ratio. Direct current amplifiers have 
been used with PbS cells with large signals but not of the order of magni- 
tude encountered in the experimental setup. Separation of the detector 
signal from the steady d.c. bias can be achieved by use of chopped radia- 
tion as indicated in Figure 5 6 An a.c. amplifier is used to amplify 
the a.c. component created by the chopper. Amplifier units were obtained 
from a Perkins-Elmer infrared spectrometer. Essentially, all the compon- 
ents except for the optical train were easily adapted for use in in in- 


frared diffraction system. The units consist of a 13 cycle per second 
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chopper, pre-anplifier, 13 cycle per second amplifier with rectified 
output, and a Leeds and Northrup recorder. With proper shielding of 
the detector, a high signal to noise ratio was obtained and relatively 


low intensity signals adequately amplified. 
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CHAPTER VI 
EXPERIMENTAL RESULTS 


Scattering data were obtained for three slices of the original 
crystal. Thicknesses were 2, 3, and 6 mm. Dislocation densities 

s e e 7 8 2 
obtained by etch-pit counts were estimated to be 10 -10 per cm for 
each specimen. Since the first tests were conducted without heat 
treating the crystals, no contribution to the light scattering by oxy- 
gen was anticipated. Microscopic examination did not reveal the 


presence of a separate phase due to SiO, clustering or other impurities. 


2 

Figures 6-1 and 6-2 represent the plotted scattering data for the 
2mm and 6 mm specimens. The 3 mm crystal produced essentially identical 
data as the 2 mm specimen, The highest value of ea correspond to 4 
scattering angle of 2 degrees. Also indicated on each curve is the limit 
specifying Rh <j 1. Theréfore the data was obtained at angles well with- 
in those for which the Guinier analysis is applicable. A check was made 
for parasitic diffraction due to the collimating slits. This was per- 
formed with the optical arrangement identical to that of the test runs 
but without the crystal. The intensity was reduced to essentially zero 
at angles of 1 degree ( aa Sharp minima were observed at an angle 
of 40 minutes. This caused the discontinuity in the specimen scattering 
curves observed at values of ae equal to 2-3. The main conclusion drawn 
from the parasitic diffraction test is the absence of any contribution to 
scattering in the region of interest, + ae due to the slits. 

From the limiting slopes, radii of gyration of 5.70 AA and 4.79 
are calculated for the 2mm and 6 mm specimens respectively. To inter- 


pret these values in terms of a particle size auxiliary information of 
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particle shape is required. As pointed out in Chapter II, the results 
of earlier investigators, suggest a poly-»zonal or circular platelet 
configuration. The difference in electron density is concentrated about 
the periphery of the loop. A good approximation of particle shape would 


be that of an annulus with a radius of gyration of: 


if? oe 
© : . a 


If the width of the dislocation line and bordering displaced atoms is 
considered to be narrow compared to the radius of the loops one can 


fesume rr, “ r For this approximation, the expression reduces to: 


i a 
In order to calculate a particle size, we must therefore assume that the 
particles are discs of negligible thickness whose radius of gyration is 
the effective radius in the plane of the disc. This assumption leads to 
a particle size with diameters of 11.4Quand 9. 28y4for the specimens of 
figures 6-1 and 6-2 respectively. These figures are of the same order 
of magnitude of the values observed by Dash using decoration techniques. 
The resistivity measurement imply a considerable impurity content. The 
presence of oxygen will be seen at large values of ae. The results in- 
dicate no particles less than 11 A, in diameter. We therefore conclude 
that the impurities are in solution in the silicon crystal or at most 


trapped in the dislocations. 





Figures 6-1] and * incicate a distribution of particle sizes. 
The long linear portion of the slope at large angles describe a larger per 
centage of smaller particles. The sharp bend in the curve at small angles 
shown in figure 6-1 indicate a smaller proportion of large particles for 
the 2 mm specimen. However, the data at angles less than 1 degree are 
not completely reliable due to the contribution of parasitic scattering 


by the collimating slits. 
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APPENDIX I 


DERIVATION OF SMALL-ANGLE SCATTERING FORMULAS 


A. Scattering By Several Waves 


A sinusoidal wave, of frequency cJ, is defined by its amplitude 
A and phase @ 
Qs ACos Cut +@} 


This can be expressed in complex variable notation: 
wi= Ae) 1fHt +8) i 
oy 
in which a is the real part of U., For several waves of same frequency 
propagated in the same direction the observed amplitude will be the sum 
of the amplitudes of each individual wave. 


J = 2 be er" sare 


The resultant amplitude and phase can be obtained by vector addition 


J i 


from the following diagram: 





Fiq I-d 


Since the intensity, I, is equal to the square of the amplitude, it is 


‘evident from the above diagram that: 


2 : 2 
I= (> A; Cos G, ) + (SA, Sin d, | 
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Expressed in complex variable nctation: 


| a ore = 
t= U0, (eset sa, e 1% \( et ZA, & -— 


CRAY An es (% ~Os) —- 22, AAK Cos (6, -@,) I-J 


B. Phase Difference of Waves Scattered by Two Electrons 








Fig BZ 


— 
In the above figure, Sy represents a unit vector along the incident rays. 


Points O and A represent two electrons scattering the incident rays. S&S 
is the scattered ray unit vector. The phase difference between the rays 


scattered at O and A can be represented by the retardation distance MA + An. 


From vector algebra: 


= 
ImAl = OA*S, 
ae ~y 
lAw| = OA a 5 - 
lmA + An| =) @.*(5,-5 


The phase difference is therefore: 


21 =—> eta 
d = N OA § ( Sos 
—_) we 
Since iS. - S{is equal to 2@Sin 9 and is parallel to ON ) 


Q = = QASin @ Cos aA ( T-2) 


4] 





Ci Intensity of Low-Anvle Scattering by Small Particles 





1. Oriented particles. The intensity of the scattering by an 
individual atom is rf where I. is the intensity of the scattering by one 


electron and f£ the atomic scattering factor. The amplitude is then re- 


presented by A= Y Te { and formula (I-1)can be written: 


1 Pesgeee Gp Za e f.. Cos &~ B,) 


D 





Fig. T-3 < 

14 O 5 
In figure I-3, € represents the scattering angle and is equal to 20. 
Since the scattering considered here is at very small angles, OD may be 


—> 
regarded as perpendicular to S for all values of € . The length ay 


is the projection of Ay on ON. ‘Therefore the phase difference, formula 


(1-2), can be written: Oo 
277 ~ aA es St A 
Q,. =~ aw € OA, Cos _~_ =.” E Oa, 





and 
= 0. eee, 
T= Te DEL, Cos( E+ Qa). 
ARE y) 
Letting Q,a. = ne and Pee Es h , the intensity can be expressed 
as a cosine series: i 4 


ae 23 E cc 
t/t, al? ee PBZ CG Fa Li) vr 24 


The various scattering points, a a ee can be projected on the 


|| 


normal OD as X 4 ee ne Hence 


aia 

i 
oN A hk geo Xx I< a x i; 
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The d b Suma tl ' " 2 : pie fj “a .- pee 
r ouble sut 101 * S : 7 fe ( Xt) Suiimed first over k is: 


Ber peo yi oe 2: x - 2 CMG, + xP 2 6, 
Choose the origin such that $ CL a » ‘This represents 
the projection of the center of gravity of the particle on OD. The 
second summation gives: 

RPe (hig) = 42,4) ° ee Sm Xk 
= A t is the total number of electrons, n, in the particle, and 
the length R, is defined such that: " 
z, s wR* = nR- * 2 oe X 
Since f, is approximately proportional to the mass of the atom, 


k 
2. 2. 
SmoeR = MR = Z Me Xx 


2 
where M is the total mass of the particle and MR, the moment of inertia 


of the particle with respect to the plane normal to OD and passing 


through the center of gravity. The double summation is then written: 


ae Fel Ri) = An R, 


and ‘ mn 
Tr, se ntL in h Ree . 5 | 


where 
ZiT 
h= 1 € 


ee : : : x 
This is the expansion of the exponential series e , Thus the formula 


for the scattering of an isolated particle is: 


eee 
r _ 2 ~ ht pe Ste i, ES 
Ite =—-~ ne 


ad Seta 
=ne Aa 
which is formula (3-3) in Chapter IIT. 
2. Random oriented particles - For randomly oriented and irregular- 


ly placed particles, the total intensity is the sum of the intensities 


diffracted by the individual particles. When there are N particles: 
Lu 


N 
. 2 > R 
a mnLi-h = 
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me 


Ea 
N 


Rye The mean value of the moment of inertia of a particle in relation 


to planes of random orientation passing through its center of gravity 





where represents the mean of the squares of the values of 


must be found. From classical mechanics: 

MRE = A (mR) +O (MRG) + POO MRE 

2 
where MR etc, are the moments of inertia with respect to the three 
principle planes and a, @ ¥* the direction cosines to the normal 


2 
of the plane considered, The mean value of Ry iG: 


MR" = = Mae at + Ry 23) = EMR* 


e. ee 
MR being the moment of inertia in relation to the center of gravity, 


R is the radius of gyration. Therefore: 


tp e 
Wr enin[i- a Pe POS | 


which is the series expansion on Rake 4n% rter 
L 2 a a mad 2 NX = » At 


This is formula(3-4) of Chapter IIT. 
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